Levy flights in confining potentials 
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We analyze confining mechanisms for Levy flights. When they evolve in suitable external po- 
tentials their variance may exist and show signatures of a superdifi'usive transport. Two classes of 
stochastic jump - type processes are considered: those driven by Langevin equation with Levy noise 
and those, named by us topological Levy processes (occurring in systems with topological complexity 
like folded polymers or complex networks and generically in inhomogeneous media) , whose Langevin 
representation is unknown and possibly nonexistent. Our major finding is that both above classes 
of processes stay in affinity and may share common stationary (eventually asymptotic) probability 
density, even if their detailed dynamical behavior look different. That generalizes and offers new 
solutions to a reverse engineering (e.g. targeted stochasticity) problem due to L Eliazar and J. 
Klafter [J. Stat. Phys. Ill, 739, (2003)]: design a Levy process whose target pdf equals a priori 
preselected one. Our observations extend to a broad class of Levy noise driven processes, like e.g. 
superdiffusion on folded polymers, geophysical flows and even climatic changes. 
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I. INTRODUCTION 

The study of random walks in complex structures is a 
key point to understanding of properties of many physi- 
cal and non-physical systems, ranging from transport in 
disordered media [1,] to transfer phenomena in biological 
cells and various real- world networks 0, Q. It is well- 
known that a mean square displacement of a freely dif- 
fusing particle depends on time linearly < X'^{t) >(x t. 
If a diffusion is anomalous, then < X'^{t) >oc t"^ , where 
7^1,0<7<oo. If7<l, the dynamics is called subd- 
ijfusive otherwise superdijfusive. A superdiffusive motion 
of a particle may be generated by means of non-Gaussian 
jump- type processes. 

At this point one often invokes Levy flights. Their 
free version may seem somewhat exotic since their sec- 
ond moments are nonexistent. However, Levy flights in 
confining external potentials show up less exotic behav- 
ior and do admit the existence of first few moments (see, 
e.g., Ref. [4|). Thus they may be employed to model a 
superdiffusive transport. 

Levy flights, being non-Gaussian jump-type processes, 
quite apart from serious technical difficulties and a short- 
age of analytically tractable examples, occur in many 
fields of modern statistical physics and have won ma- 
jor attention in the last two decades [3|-[1H- Most of the 
current research is devoted to Langevin equation based 
derivations, where a deterministic force is perturbed by 
the (white) noise of interest, [T^-[23|. However, in a 
number of publications, another class of jump-type pro- 
cesses was introduced under the name of topologically 



'Electronic address: pgarQuni. opole.^ 
t Electronic address: jStef@um. opole. pl| 



induced super-diffusions, The origin of this name 

is due to the fact that such processes occur primarily 
in the systems with topological complexity like folded 
polymers or complex networks. An observation of 0] 
was that topological super-diffusion processes do not por- 
tray a situation equivalent to any of standard fractional 
Fokker-Planck equations and seem not to correspond to 
any Langevin equation. On the other hand, in the discus- 
sion of above topological Levy processes main emphasis 
has been put on their super-diffusive behavior with some 
neglect of confining effects and the resultant emergence 
of non-Gibbsian stationary probability densities, [sj-f^. 

We address the latter issue and set general confine- 
ment criteria for an analytically tractable case of Cauchy 
noise-driven processes. The results obtained appear to 
be more general and not specific to Cauchy noise. To 
this end, some ideas have been adopted from the gen- 
eral theory of diffusion-type stochastic processes where 
an asymptotic approach towards equilibrium (stationary 
probability density function (pdf)) is one of major topics 
of interest, [221] . 

To handle topological Levy processes we use a conve- 
nient and general mathematical tool, named Schrodinger 
(or Levy - Schrodinger for non-Gaussian processes) semi- 
group. This tool naturally appears if one attempts to 
transform the evolution equation for the pdf p of a certain 
stochastic process (e.g. standard or fractional Fokker- 
Planck equation), into the time - dependent Schrodinger- 
type equation (the parabolic one in the Gaussian con- 
text; there is no imaginary unit before time derivative) 
dttp — Ti-tp- Here, Ti. receives a natural interpretation 
of a Hamiltonian operator, — 7i stands for a semigroup 
generator. A proper exploitation of a semigroup opera- 
tor exp(— tTi) allows not only to generate the evolution 
equation for the pdf (differential or pseudo-differential in 
case of non-Gaussian Levy noise, see below) but gives 
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access to hitherto unexploited evolution scenarios which 
are not captured by the standard Langevin modehng. 

We shall demonstrate that topologically induced pro- 
cesses of Refs. 0-0 form a subclass of its solutions 
with a properly tailored dynamical semigroup and its 
(Feynman-Kac) potential, That allows to take 

advantage of the existing mathematical theory of Levy 
processes and Levy - Schrodinger semigroups, [3, 
and [1, [m, , where free Levy noise generators are ad- 
ditively perturbed by suitable confining potentials. The 
theory works well for both Gaussian and non-Gaussian 
processes. 

We note here, that in the Brownian case, the 
Schrodinger problem incorporates the well known trans- 
formation of a Fokker-Planck equation into a generalized 
diffusion equation, e.g. a transition to the Hermi- 
tian (strictly speaking, self-adjoint) problem whose eigen- 
function expansions yield transition pdfs of the pertinent 
process. 

In this article, we consider an impact of external confin- 
ing potentials upon Levy flights. The flights may be influ- 
enced directly or indirectly (here via conservative forces) 
leading to inequivalent Levy processes. An indirect influ- 
ence refer to Langevin modeling, while a direct one refers 
the Levy semigroups. While making this specific distinc- 
tion between the two ways of response of Levy noise to 
external potentials, we address an issue of an apparent 
incompatibility between them, raised earlier [3] ■ The re- 
sults obtained set a bridge between these seemingly dif- 
ferent classes and may shed some light on the emergence 
of varied types of a superdiffusive dynamics in complex 
structures, especially those involving significant spatial 
inhomogeneities. 



II. THEORETICAL FRAMEWORK 

A. Smoluchowski processes and Schrodinger 
semigroups 

To make paper self-contained, here we recapitulate 
the main derivations, which will be necessary for us in 
subsequent discussion. We begin with consideration of 
a one- dimensional (ID) Smoluchowski diffusion process 
p^ . with the Langevin representation x = b{x, t) + A{t), 
where {A{s)) = 0, {A{s)A{s')) = 2D5{s - s'). Here, 
b{x, t) is a forward drift of the process, admitted to be 
time - dependent, unless we ultimately pass to Smolu- 
chowski diffusion processes where b{x, t) = b{x) for all 
times. 

If an initial pdf pq{x) is given, then the diffusion pro- 
cess drives it in accordance with the Fokker-Planck equa- 
tion dtp = DAp — V {bp) (in the ID case V = d/dx, 
A = d'^/dx^). We introduce an osmotic velocity field 
u = 13 In p, together with the current velocity field 
V = b — u. The latter obeys the continuity equation 
dtp = — Vj, where j — v ■ p has a standard interpretation 
of a probability current. The time-independent drifts 



b{x) of the diffusion processes are induced by external 
(conservative, Newtonian) force fields / = —W. One 
arrives at Smoluchowski diffusion processes by setting 



J_ 

to/3 



1 

to/3 



VV. 



(1) 



Here, to is a mass and /3 is a reciprocal relaxation time 
of a system. The expression ([T]) accounts for a fully - 
fledged phase - space derivation of the spatial process, 
in the regime of large (3. It is taken for granted that 
the fluctuation-dissipation balance gives rise to the stan- 
dard form D = kBT/mp of the diffusion coefficient D (T 
stands for a temperature and ks is Boltzmann constant). 

Let us consider a stationary asymptotic regime, where 
7 J* = 0. We denote an (a priori assumed to exist, 
[22|), invariant pdf = p*{x). Since in stationary case 
V — = 0, we have 



L»Vlnp*. 



(2) 



Since b = f /mP does not depend on pdf explicitly, 6 = 
6* and p^{x) — (l/Z) exp[—V{x)/kBT]. It is seen that 
our outcome has Gibbs-Boltzmann form with Z being a 
partition function, Z = J cxp(— y/fc^T) dx. 
Denoting = —kBTln Z , we have 

p4x) - exp ([F, - V{x)]/kBT) = exp[2$(x)] . (3) 

Here, to comply with the notations of Refs. [3? ] and 
with subsequent discussion of a topological generalization 
of the Brownian motion and then Levy flights we 
have introduced a new potential function such that 



P 



^1'^ - exp($) and b = 2DV'^. 



Following a standard procedure [ld | we transform 
the Fokker-Planck equation into an associated Her- 
mitian problem by means of redefinition p(x, t) — 
6* {x,t) exp[^{x)], that takes the Fokker-Plack equation 
into a parabolic one 9*6** — — Its po- 

tential V derives from a compatibility condition V{x) = 
{l/2W/{2D) + Vb]. 

Smoluchowski process with a unique asymptotic Gibb- 
sian pdf implies 



V = D 



1/2 



1/2 



(4) 



This equation is a trivialized version (due to the time- 
independence of its solution) of the time adjoint equation 

dtO = -DM + V0, see Refs. H [HI setting 9 = p^^ 

Introducing {\/2mD rescaled) Schrodinger- type 
Hamiltonian Ti. — —DA. + V, one arrives at a dynamical 
(Schrodinger) semigroup operator exp(— iTi), with the 
dynamical rule 6*{t) — [exp(-i7Y)(?*](0), taking forward 
the initial data 6*[x, 0). 

For completeness of discussion, we note that the time 
adjoint equation, if applicable, would come out from 
the reverse time evolution taking a given final (ter- 
minal) 0{x,tfin) backwards in time to 6{x,tfin ~ t) = 
[exp(— i7i)0](ifin), all motions being confined to an in- 
terval [0, ifin]. 
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B. Levy - Schrodinger semigroups. 

Before passing to an analysis of Levy flights, let us set 
general rules of the game with respect to the response to 
external potentials, once a free noise is chosen. We re- 
call that a characteristic function of a random variable X 
completely determines a probability distribution of that 
variable. If this distribution admits a pdf p{x), we can 
write < exp(zpX) >= j ^ p{x) exp{ipx)dx which, for in- 
finitely divisible probability laws, gives rise to the famous 
Levy-Khintchine formula (see, e.g. [l3|) 



< exjp{ipX) >~ cxp< iap — (cr^/2)p^ 



refers to the Cauchy process, see e.g. d, [TT|, [13] ■ The 
pseudo - differential Fokker-Planck equation, which cor- 
responds to the fractional Hamiltonian (28) and the frac- 
tional semigroup exp{—tH^) = cxp(— A| A|'^/^), reads 



(9) 



to be compared with the conventional heat equation 
dtp = DAp. 

For a pseudo-differential operator |A|''/^, the action 
on a function from its domain is greatly simplified (as 
compared to Levy-Khintchine formula ([HI)), in view of 
the properties of the Levy measure v^j^{dx). We have 

B a [HE Hi: 



exp(ipy) - 1 - 



v{dy)), (5) 



where i'{dy) stands for so-called Levy measure. By 
disregarding the deterministic and jump-type contribu- 
tions in the above, we are left with < exp(zpa;) >= 
exp(-(TV/2), hence p{x) = exp{-x^ /2a'^). 

In terms of the random variable Xt ~ {2Dy/^At of the 
Wiener process, we have < exp{ipXt) >— exp{~tDp'^). 
By employing p ^ p = — iV we identify the semigroup 
operator exp{tDA), with A = d^ /dx"^. This involves a 



Dp^ 



-DA of the general Hamil- 



special version Ti 
tonian H = F{p). 

From now on, we concentrate on the integral part of 
the Levy-Khintchine formula, which is responsible for ar- 
bitrary stochastic jump features. By disregarding the 
deterministic and Brownian motion entries we arrive at: 



F{p) = 



exp{ipy) — 1 — 



ipy 



1 + r 



where v{dy) stands for the appropriate Levy measure. 
The corresponding non-Gaussian Markov process is char- 
acterized by < exY>{ipXt) >— exp[— ti^(p)] and yields an 
operator F{p) = H, with p = —i\7. 

For the sake of clarity we restrict further considera- 
tions to non-Gaussian random variables whose pdf 's are 
centered and symmetric, e.g. a subclass of stable distri- 
butions characterized by 



i^(p) = A|p|^ = A|A|''/2_ 



Here fi < 2 and A > stands for the intensity parame- 
ter of the Levy process. The fractional Hamiltonian H, 
which is a non-local pseudo-differential operator, by con- 
struction is positive and self-adjoint on a properly tai- 
lored domain. A sufficient and necessary condition for 
both these properties to hold true is that the pdf of the 
Levy process is symmetric, |14|. 

The associated jump-type dynamics is interpreted in 
terms of Levy flights. In particular 



F{p) = Abl ^n = F{p) = A|V| = A(-A) 



1/2 



{\A\^'^f){x) 



[f{x + y)- f{x)]y^{dy). (10) 



The Cauchy-Levy measure, associated with the Cauchy 
semigroup generator |A|^/^ = |V|, reads 



vi/2[dy) = 2 ■ 



(11) 



The substitution y ^ z = x + y permits to reduce the 
Eq. pn)) to the familiar form 



(|V|/)(x) 



1 



/(^) - /(^) 
\z — 



dz, 



(12) 



where l/7r|z — has an interpretation of an intensity 
with which jumps of the size \z — x\ occur. 



III. RESPONSE TO EXTERNAL POTENTIALS: 
^{dy), (6) STATIONARY DENSITIES 



A. Langevin modeling 

The pseudo-differential Fokker-Planck equation, which 
corresponds to the fractional Hamiltonian ([7]) and the 
fractional semigroup exp(— iTi^) = exp(— tA| A|^/^), has 
the form ^ to be compared with the Fokker-Planck 
equation for freely diffusing particle (or above heat trans- 
fer equation) dtp = DAp. 

In case of jump- type (Levy) processes a response to ex- 

(7) ternal perturbations by conservative force fields appears 
to be particularly interesting. On one hand, one encoun- 
ters a widely accepted reasoning (Refs. [131- US]) where 
the Langevin equation, with additive deterministic and 
Levy "white noise" terms, is found to imply a fractional 
Fokker-Planck equation, whose form faithfully parallels 
the Brownian version, e.g. (c.f. Ref. [13, see also ^) 

X = h{x) + A^'{t) =^ dtp ^ -V(6 • p) - A| Aj^/^p . (13) 

Here we make a remark regarding our notations. In ID 
case operator V means simply differentiation over x (see 

(8) also above) so that all quantities like / are scalars. In 
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higher dimensions the operator V, acting on vector quan- 
tity b-p {b = — Vy/m/3) should be understood as a vector 
divergence, i.e. the term Vib ■ p) = div(6 • p). Also, here 
we emphasize a difference in sign in the second term of 
Eq. dlS]) as compared to that in Eq. (4) of Ref. [13] ■ 
There, the minus sign is absorbed in the adopted defi- 
nition of the (Riesz) fractional derivative. Apart from 
the formal resemblance of operator symbols, we do not 
directly employ fractional derivatives in our formalism. 

B. Topological route 

The other approach to account for external perturba- 
tions is that, by mimicking the above Gaussian strategy, 
we can directly refer to the Hamiltonian framework and 
dynamical semigroups with Levy generators being addi- 
tively perturbed by a suitable potential. For example, as- 
suming that the functional form of V{x) guarantees that 
Tifj^ = A|A|''/^ -I- V is self-adjoint and bounded from be- 
low, we may pass to the fractional (non-Gaussian, jump 
process) analog of the generalized diffusion equation: 

dte, = -x\/:^\t'''^e,-ve,. (i4) 

The dynamical semigroup reads exp(— tTi^) and the com- 
patibility condition related to Eq. ([3]), takes the form of 
the time-adjoint equation dtO = A| Ap'/^gi^ V6I Gen- 
eral theory [ill tells us that 9* [x,t)9{x,t) = p(x,t) 



We note a systematic sign difference between our |A|'^/^ 
and the corresponding fractional derivative A^/^ of Refs. 

[i,i[i. 



C. A discord and the reverse engineering problem 

The puzzling point is that for the Levy process in ex- 
ternal force fields, the Langevin approach yields a conti- 
nuity (e.g. fractional Fokker- Planck) equation in a very 
different form 

9*P=-vf-^p) -A|Ar/2p. (18) 
\ mp J 

The conclusion of Refs. [1|-[1] was that, while assuming 
^ V where V is (up to inessential factors) the above 
external force potential, the two transport equations Ijl6p 
and are plainly incompatible so that Eq. seems 
not correspond to any Langevin equation with Levy noise 
term and b = —W/m(3 as a deterministic part and vice 



stands for a pdf of an affiliated Markov process that inter- 
polates between the boundary data p{x, 0) and p(x, ifin)j 
at times t £ [0, ign]. 

We consider time- independent 9{x,t) = 6{x) and 
hereby mimic the Gaussian ansatz: 9{x) = exp[$(a;)] so 
that 9*{x,t) = p{x,t) exp[-$(x)]. If we set exp[<l5(a;)] = 

1 /2 

pj (x), we get the compatibility condition (see Eq. Q): 

V = -a I^I f* . (15) 
p* 

This identity should be compared with Eq. (8) in Ref. 
^8], where an analogous effective potential was deduced 
in the study of Levy flights in inhomogeneous media. 

In view of the semigroup dynamics, we deduce a con- 
tinuity equation with an explicit fractional input 

dtp = 9dt9* = -A(exp$)|A|^/2[gxp(-$)p] +v-p. (16) 

Up to cosmetic changes $ —V/2kBT (compare with 
Eq. ([3])), Eq. is identical with transport equations 
employed in a number of papers. There, the investigated 
process was named a topologically induced superdiffu- 
sion. Namely, with respect to explicit form of Eq. (|15p . 
the Eq. (|16p assumes a familiar form of the transport 
equation (with respect to A = 1 and k = l/iksT)), see 
Eq. {&) in Ref. [1], Eq. (5) in Ref. [| and Eq. (36) in 
Ref. H. 



(17) 

I 

versa. This puzzling discrepancy has not been explored 
previously in more depth. 
The problem we address is: 

(i) choose a functional form of V{x) and thus the drift 

of the Langevin - type process; 

(ii) infer an invariant pdf that is compatible with the 

fractional Fokker-Planck equation ; 

(iii) given p*, deduce the Feynman-Kac (e.g. dynamical 
semigroup) potential V by means of Eq. (jlSp ; 

(iv) use V in p6)) and verify whether the "topologi- 
cally induced dynamics" is at all related to that 
associated with p^ (and thus to the underlying 
Langevin equation with Levy noise); 

(v) check an asymptotic behavior of p(x, t) in both sce- 

narios p6)) and pS]) to find possible differences in 
the speed (convergence time rate) with which the 
common invariant pdf p^{x) from item (ii) is ap- 
proached; 



I 

dtp = -exp(-KF/2) |A|'^/2exp(KT//2)p + pexp(KP"/2)|A|''/2 exp(-KF/2), 
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(vi) repeat the procedure in reverse order by starting 
from step (iii) and then deduce the drift for the 
Langevin equation with Levy noise; next compare 
the dynamical scenarios (fTB)) and for any com- 
mon initial pdf. 

We recall that the the above problem is non-existent in 
the case of Brownian motion. There, the Fokker-Planck 
dynamics and the related parabolic equations do refer to 
the same diffusion-type process. 

We shall demonstrate below that both Langevin - 
driven and semigroup - driven Cauchy processes, albeit 
non - coinciding literally, keep resemblance to each other 
and may share common for both stationary pdf. A su- 
perdiffusive dynamical behavior is generically expected 
to arise and an asymptotic approach towards a station- 
ary pdf is then in principle possible. This motivates the 
" targeted stochasticity" discussion whose original formu- 
lation (in terms of the reverse engineering problem) for 
Langevin - driven Levy systems can be found in Ref. [2^ . 
The original formulation of the reverse engineering prob- 
lem reads: given a stationary pdf, can we tailor a drift 
function so that the system Langevin dynamics would 
admit the predefined as an asymptotic target? 

We employ the reverse engineering problem to analyze 
Cauchy processes in confining potentials. In the course of 
the discussion, we in fact extend its range of applicability 
(that applies to more general stable processes as well) 
and demonstrate that a priori chosen stationary pdf may 
serve as a target density for both Langevin and semigroup 

- driven Cauchy processes. Even though their detailed 
dynamical patterns of behavior are different. In the near 

- equilibrium regime this dynamical distinction becomes 
immaterial. 



Fokker-Planck equation dtp = — A| V|p -I- V[(7a;)/9] reads 



P*{x) 



X 

7 



(19) 



c.f. Eq. (9) in Ref. [ij]- Here, the modified noise 
intensity parameter a is a ratio of an intensity param- 
eter A of the Cauchy noise and of the friction coeffi- 
cient 7. Note that a characteristic function of this pdf 
reads F(p) = —cf\p\ and accounts for a non-thermal 
fluctuation-dissipation balance. 

For Cauchy random variable Xt we have {exp{ipXt)) 
= exp(tA|p|). The corresponding pdf has the form (fT5)) 
with a ^ tX, e.g. p{x,t) — Xt/TT[{Xt)'^ + x'^]. Here, a and 
tX play a role of scaling parameters specifying the half- 
width of the Cauchy pdf at its half-maximum. Since tX 
grows monotonically, the free Cauchy noise pdf flattens 
and its maximum drops down in time. 

Since a — A/7, the confining drift —7a; may stop the 
"flattening" of the probability distribution and stabilize 
the pdf at quite arbitrary shape (with respect to its max- 
imum and half- width, see above), by manipulating 7. For 
example, 7^1 implies a significant shrinking of the dis- 
tribution as compared to the reference (free noise) pdf 
at any time t ^ 1/A. In parallel, a maximum pdf value 
would increase: I/ttA — > 7/7rA. 

The OUC case refers to Cauchy flights in a confining 
(harmonic) potential, but does not imply the confined 
flight, since the variance of the asymptotic density di- 
verges. We note that confined Levy flights and specifi- 
cally confined Cauchy flights, have been analyzed earlier 
in Refs. [H-ll^. 

To deduce the potential V for the OUC process with 
given invariant pdf p*, we need to evaluate the right- 
hand-side of the defining Eq. (|15p . with p = 1. We 
employ Eq. (|12p . so arriving at: 



IV. CAUCHY DRIVER 

In view of serious technical difRculties we shall not at- 
tempt to present a fully fledged solution to the above 
formulated problem for any symmetric stable jump-type 
process and any conceivable drift. Instead, we turn our 
attention to situations where explicit functional forms of 
invariant densities are available. Most of them were in- 
ferred in the problems, related to Cauchy noise, see Refs 
[i3"[23]- In particular, attention has been paid 
to confining properties of various drifts upon the Cauchy 
noise. On the other hand. Levy flights through a "po- 
tential landscape" (topological processes of Refs. [1]-@|) 
were interpreted as (enhanced) super-diffusions. 



A. Ornstein - Uhlenbeck - Cauchy process 

Let us consider the Ornstein - Uhlenbeck - Cauchy 
(OUC) process, whose drift is given by b{x) = —^x, and 
an asymptotic invariant pdf associated with the Cauchy- 



1 



A ((72 +a;2)i/; 



■V[x) 



^cr2 + (x + y)2 ^/ a'^ + : 



%■ (20) 



Because of the integrand singularity at ?/ = 0, we must 
handle the integral in terms of its principal value. Intro- 
ducing the notation a ~ + , arrive at, [2^ : 



V(.x) = - 

TT 



2 X Jli + X 

^ + - In ^ 

^ a a — x 



(21) 



Here, V(a;) is bounded both from below and above, with 
the asymptotics (2/|a;|) In at infinities, well fitting to 
the general mathematical construction of (topological) 
Cauchy processes in external potentials, see Ref. hj for 
details. The plot of potential (pij) is reported in Fig. 1 . 

Accordingly, we know for sure that there exists a topo- 
logical Cauchy process with the Feyman - Kac potential 
V, Eq. (HH), whose invariant density coincides with that 
for the Langevin - supported OUC process. 
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FIG. 1: The coordinate dependence of potentials V(a;): (|2ip for different a (main left panel), (|24|) (inset to left panel) and (|35|) 
for different /? (right panel). 



B. Confined Cauchy processes: Langevin and 
topological targeting 

To analyze a time-dependent behavior of both topo- 
logical and Langevin - driven process, below we con- 
sider specific numerical example, admitting finite vari- 
ance < X'^it) >. This time dependent variance permits 
to analyze a particular scenario of approaching the in- 
variant (equilibrium) density in the large time regime. 
We will see, that two considered jump - type processes, 
whose time evolution is embodied respectively in the frac- 
tional Fokker- Planck equation and in Levy-Schrodinger 
semigroup (topological case) dynamics are definitely alike 
as they share a common invariant density. In the near- 
equilibrium regime, any dynamical distinction between 
these motion scenarios becomes immaterial. However, 
their detailed dynamical behavior far from equilibrium 
might be different and this issue deserves further analyt- 
ical and numerical exploration. 

To our current knowledge, there is no Langevin - type 
representation of a topological process and vice versa, 
even though an invariant density is common for both. 
Nonetheless, we will demonstrate that by starting from 
a common initial probability density, the two (Langevin 
and dynamical semigroup) motion scenarios end up at a 
common invariant density. 

Neither OUC process nor its topological counterpart 
are confined. For the Cauchy density, the second mo- 
ment is nonexistent. We shall verify the outcome of the 
OUC discussion for Cauchy-type processes whose invari- 
ant densities admit the second moment due to confine- 
ment. Let us consider the quadratic Cauchy pdf: 



1 



TT (1 + X^) 



2\2 



(22) 



Now, let us proceed in reverse order departing from Eq. 
so that {l/V^)pV^ = (l/7r)/(l -t- a;^) is actually 



1/2 

Cauchy pdf. We consider f{x) = pj as the initial data 
for the free Cauchy evolution dtf — A|V|/. This takes 
f{x) into the form 



l + Xt 



TT {l + Xt)-^+X^ 



Since A| V|/ = — limj^o dtf we end up with 
limt^oft/ a;2 - 1 

v(x) = — J — = 



(23) 



(24) 



The shape of this potential is shown in Fig.l (inset to 
left panel). A minimum —A is achieved at a; = 0, V = 
occurs for a; = ±1, a maximum +A is reached at x — > 
±oo. 

The potential is bounded both from below and above 
and hence can trivially be made non-negative (add A). 
This means that the potential (|24p is fully compatible 
with the general construction of Ref, [llj. This topo- 
logical process is generated by Cauchy generator plus a 
potential function, see Ref. [ll|, is of the jump-type and 
can be obtained as an e — > limit of a step process with 
a minimal step size e. 

Note, that in Ref. [11] no explicit example of the con- 
fining potential V has been proposed. Eqs. (PT|) . ([M]) 
provide such examples, which, to our knowledge, have 
never been exploited in the literature. 

At this point, let us make a guess that the quadratic 
Cauchy pdf actually stands for an invariant pdf of the 
" normal" Langevin-based fractional Fokker-Planck equa- 
tion ^TE\i with a drift of the form ([T|). Accordingly we 
should have dtp* = = — V(6p*) — 7|V|p* and therefore 
the admissible drift function, if any, may be deduced by 
means of an indefinite integral: 



bix) 
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{\V\p,){x)dx. 



(25) 
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For quadratic Cauchy pdf (I22p the explicit form of b{x) 
H reads 



6(x) = -^(:e2^3). 



(26) 



Thus, there exists the Langevin process whose invariant 
pdf is shared with a corresponding topological process. 
In the near - equilibrium regime a dynamical distinction 
between the pertinent processes becomes immaterial. In 
other words, if we wish to deal with the Langevin process 
associated with the quadratic Cauchy density (|^^ . the 
proper drift form is given in (j26p . 

To analyze numerically the above apparent discord be- 
tween Langevin-driven and topological processes, we use 
the invariant pdf ((22l) . having drift (l26l) and Feynman - 
Kac potential ([M)) . We have chosen this invariant pdf 
as it has a finite variance, which permits us to capture 
the details of near - equilibrium, initial and intermediate 
stages of time evolution. 

For numerical solution we use simple Euler scheme for 
time derivatives and numerical integration (more specif- 
ically, we calculate Cauchy principal value of integrals) 
on the each Euler time step for evaluation of fractional 
derivative |V|. The initial state corresponds to a parti- 
cle localized at a; = 0, corresponding to the minima of 
both potential, derived from the drift (|26p and Feynman 
- Kac potential p{x,t = 0) = 6{x). The solutions 

p{x,t) of the equations (Langevin-type process) and 
p?)) (topological process) are reported in Fig. 2 (upper 
and middle panels respectively). 

It is seen that topological diffusion process needs more 
time to achieve the invariant pdf, appears to be slowed 
down as compared to the Langevin scenario. This is il- 
lustrated in Fig. 2 (lower panel), where the time evolu- 
tion of variances for both processes have been plotted. 
The time evolution occurs from zero variance of 5 func- 
tion to asymptotic variance < X'^{t oo) >= 1 of the 
pdf (P^ . It is seen, that variance for Langevin - type 
process achieves the asymptotic value at (dimensionless) 
time t « 0.5, while for topological diffusion this time 
t « 2. The shapes of < X'^{t) > for both processes defi- 
nitely resemble a super - diffusive motion. 



C. Confined Cauchy family 
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Now we consider a broader class of pdf 's related to the 
Cauchy noise. Any continuous pdf p can be associated 
with Shannon entropy S{p) = — J p In pdx, t25j. If an 
expectation value < ln(l + x'^) > is fixed, the maximum 
entropy probability function belongs to a one-parameter 
family 

^* ^""^ = v^nl-1/2)) (i+x2)- (^^^ 



FIG. 2: Time evolution of pdf's p{x,t) for topological (up- 
per panel, A = 1) and Langevin-type (middle panel, [3 = 1, 
m = 1) processes. The common equilibrium pdf p* (|22p is 
also shown. Lower panel reports the time dependent variance 
X^it) for Langevin-type (solid line) and topological (dashed 
line) processes. Points correspond to numerical calculation, 
lines are guides for eye. 



where a > 1/2, [H]. 
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Cauchy distribution is a special case of the above 
that corresponds to a = 1. The density ([^^ is the sec- 
ond, a ~ 2, member of the a - integer hierarchy (we 
assume that a = 1). 

Our tentative analysis shows that for integer and half - 
integer a, the invariant pdf (f27|) admits V{x), which fits 
the restrictions of Corollary 2 in Ref. [ll| . The question 
about arbitrary a is still under investigation. 

For each specific function V{x), the resulting Markov 
jump-type stochastic process, determined by the Cauchy 
generator plus a suitable potential function, appears to 
be unique. Here we present only one specific example, 
namely we consider 



16 

P*(x) = — 



1 



57r (1 + a;2)4 



(28) 



Substitution of Eq. ([55]) into Eq. p5|) with respect to 
definition l|12p yields [2J| the following expression for the 
Feynman-Kac (semigroup) potential 



V(x) = 



J X* + 6x^ 



1 



(29) 



The potential is bounded from below, its minimum at 
X = equals —87/2. For large values of |a;|, the poten- 
tial behaves as ^ {j/2)x'^ i.e. demonstrates a harmonic 
behavior. 

Apart from the unboundedness of V{x) from above, 
this potential obeys the minimal requirements of Corol- 
lary 2 in Ref. [ll|: can be made positive (add a suit- 
able constant), is locally bounded (e.g. is bounded on 
each compact set) and is measurable (e.g. can be ap- 
proximated with arbitrary precision by step functions se- 
quences). The Cauchy generator plus the potential pS]) 
determine uniquely an associated Markov process of the 
jump-type and its step process approximations. 



Having the density ((28|) . we can readily address the 
problem (vi) of Section llll CI Namely, inserting Eq. 
to Eq. ([25| . we obtain 



b{x) 



72: 



{5x^ + 21x'^ + 35x^ + 35). 



(30) 



This function shows a linear friction b ^ —x for small x 
and a strong taming behavior b ^ —x"^ for large x. 
Let us finally consider a bimodal pdf (see, e.g. Ref. 

M) 



P*{x) 



(3^ 



1 



TT x^ - (3^X^ + (3'^' 



(31) 



which is a solution of so-called quartic Cauchy oscilla- 
tor. As a form of the (confining) potential V{x) oc x** is 
known for that pdf, we can check the correctness of the 
procedure ([55]) of deriving a drift (and hence the potential 
V{x) in Langevin scenario) for this pdf. The application 
of operator (fT^ to function (PT]) yields 



|V|p,(a;) 



a;4 _^ p2^2 _ 3^4 



(32) 



which after integration over x and division over p*(a; 
([3T|) yields 



V{x) = - / b{x)dx = -^x'^, 



(33) 
(34) 



which is exactly the form of the potential for quartic 
Cauchy oscillator. The expression pT|) can also be used 
to calculate the "topological" potential V{x) 



TT 



1 

/3^x^ + 13^ 



(35) 



Since an analytic outcome has proved not to be tractable, 
we have reiterated to numerics. The result of numerical 
calculation of the function ([35| is reported in Fig 1 (right 
panel) for different (3. ft is seen that this potential is 
also bounded from below and above, can be made non - 
negative and have all properties imposed by Corollary 2 
of Ref. [ll|- 



V. CONCLUSIONS 

Explicitly solvable models are scarce in theoretical 
studies of Levy flights, especially in the presence of 
external potentials and/or external conservative forces. 
Therefore, our major task was to find novel analytically 
tractable examples, that would shed some light on ap- 
parent discrepancies between dynamical patterns of be- 
havior associated with two different fractional transport 
equations that are met in the literature on Levy flights. 

Although the predominant part of this research is 
devoted to the standard Langevin modeling, we have 
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demonstrated that so - called topological Levy processes 
form a subclass of solutions to the Schrodinger boundary 
data problem. The pertinent dynamical behavior stems 
form a suitable Levy - Schrodinger semigroup. The cru- 
cial role of the involved Feynman - Kac potential has 
been identified. We have explicitly derived these poten- 
tial functions in a number of cases. 

The major gain of above observations is that a math- 
ematical theory of Ref. tells one what are the nec- 
essary functional properties of admissible Feynman-Kac 
potentials. Their proper choice makes a topological Levy 
process a well behaved mathematical construction, with 
a well defined Markovian dynamics and stationary pdf. 

Our focus was upon confinement mechanisms that 
tame Levy flights to the extent that second moments of 
their probability densities exist. We have shown that the 
dynamical behavior of both above classes of processes are 
close to each other in the near-equilibrium regime and ad- 
mit common (for both classes) stationary pdf. This pdf, 
in turn, determines a functional form of the aforemen- 
tioned (semigroup defining) potential function. 

We have generalized the reverse engineering (targeted 
stochasticity) problem of Ref. [l^ beyond the original 
Levy - Langevin processes setting. We have demon- 
strated, that within the targeted stochasticity frame- 
work, the concept of Levy flights in confining potentials 
is not limited to the standard Langevin scenario. The 
Levy - Schrodinger semigroup explicitly involves confin- 
ing potentials, but with no obvious link to a Langevin 
representation. Our version of the reverse engineering 
problem amounts to reconstructing from a given (target) 
stationary density the potential functions that either: (i) 
define the forward drift of the Langevin process, or (ii) 
enter the Schrodinger - type Hamiltonian expression in 
the semigroup dynamics. Both dynamical scenarios are 
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